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Abstract. The most general form of a marginal extended perturbation in a two- 
dimensional system is deduced from scaling considerations. It includes as particular 
i cases extended perturbations decaying either from a surface, a line or a point for 

jy^ , which exact results have been previously obtained. The first— order corrections to 

the local exponents, which are functions of the amplitude of the defect, are deduced 
from a perturbation expansion of the two-point correlation functions. Assuming 
covariance under conformal transformation, the perturbed system is mapped onto a 
cylinder. Working in the Hamiltonian limit, the first-order corrections to the lowest 
gaps are calculated for the Ising model. The results confirm the validity of the gap- 
^ ' exponent relations for the perturbed system. 

ON ' 

1. Introduction 

in ■ 

Following their application to semi-infinite two-dimensional systems [1], the methods 
of conformal invariance have been used to study the local critical behaviour near 
line defects [2-4] and star-shaped defects [5] in the two-dimensional Ising model. 
Although the symmetries usually considered to be necessary for conformal invariance 
[6, 7] (dilatation, rotation and translation invariance) are partially broken by such 
perturbations, the properties associated with conformal invariance are preserved and 
' the spectrum-generating algebra have even been identified [8-10]. 

, There are some indications that scaling operators remain covariant under con- 

formal transformation in the presence of marginal extended perturbations too. This 
type of defect, which was originally introduced as a deviation from the bulk coupling 
strength with a power-law decay from a free surface in the 2D Ising model [11], has 
been extensively studied in the last years [12-17]. It becomes marginal when the 
decay exponent is equal to the scaling dimension of the coupling (see section 2 below). 
Then local exponents may vary with the perturbation amplitude. When the system 
is mapped onto a strip, provided the perturbation profile is properly transformed, the 
gap-exponent relations and the tower-like structure of the spectrum are preserved 
[18-20] but the associated algebra has not yet been determined. Similar results were 
obtained with an extended perturbation induced by an internal line defect [21]. This 
type of perturbation was first considered in [22] (see also [23-25]). In the case of a 
radial extended perturbation [26, 27] or with a decaying surface field [28], the gap- 
exponent relations are also satisfied although the equidistant-level structure of the 
spectrum is lost. 
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In the present work we consider the most general marginal extended perturbation 
in the 2D Ising model and show that, up to first order in the defect amplitude, the gap- 
exponents relations remain valid. The form of the perturbation is obtained through 
scaling considerations in section 2. The first-order corrections to the local exponents 
are deduced from a perturbation calculation of the two-point correlation functions 
to logarithmic accuracy in section 3. This is followed by a determination of the 
transformed profile in the strip geometry and a study of the condition for shape 
invariance under the special conformal transformation in section 4. The first-order 
shifts of the lowest gaps are calculated and the validity of the gap-exponent relations 
examined in section 5. Specific examples are discussed in section 6. 



2. Scaling considerations 

In a continuum description let the Hamiltonian of a two-dimensional inhomogencous 
system be written as 

-0H = -0H C + J A(r)<f>(r)d 2 r (1) 

where H c is the critical point Hamiltonian of a conformally invariant system, A(r) is 
an extended perturbation, with scaling dimension y^, which is conjugate to the local 
operator <j){r) with bulk scaling dimension = 2 — y^. This perturbation can be 
written as the product of an amplitude g by a shape function Z(r) which gives the 
form of the inhomogeneity. It is assumed to be covariant under rescaling so that, in 
polar coordinates, 

z(^,e)=trz(r,e) (2) 

where the scale invariance of the angle 8 has been taken into account. With b = r 
equation (2) leads to a power-law behaviour for the radial part 

ZM) = M (3) 

whereas the angular dependence remains arbitrary. 
The inhomogeneity A(r) transforms according to 

A'(r') = g'Z (*-,0) = b y *A(r) = V>*- u gZ (4) 

As a consequence, the perturbation amplitude scales like [29,30] 

9 = W*~"g (5) 

When lo > y^ the perturbation decays strongly enough for the extended perturbation 
to be irrelevant. When u < y^, on the contrary, the amplitude increases under 
rescaling, the perturbation is relevant and the original fixed point is unstable. More 
interesting is the borderline situation where u> — y^. Then the extended perturbation 
is marginal and non-universal ((/-dependent) local exponents are expected. These 
results are easily generalized in higher dimensions. 
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3. Perturbation theory 

We now specialize to the case where the perturbation term in (1) is marginal and 
involves the energy density operator e(r). Then 

-(3H = -0H C +g J Z(r)e(r)d 2 r (6) 

The first-order change in the local critical behaviour can be deduced from a 
perturbation expansion of the correlation functions [22]. The order parameter two- 
point correlation function has the following expansion in powers of g 

G aa (R,g) = £ 9 - / ((a(0)a(R)e(r 1 )s(r 2 ) ■ ■ ■ ^J» II Z ( r i) d \ ( 7 ) 

n=0 H - J i=l 

where the double brackets denote the irreducible part of a multi-point correlation 
function and R is far from the defect in the bulk of the system. The n-th order 
contribution can be rewritten as 

6G%>=g n / ((o-(0)*(R)e(r 1 )s(r 2 )...s(r n )))l[Z(r i )d 2 r i (8) 

Jri<ro---<r„ „-_i 



i=l 



where the n! in the prefactor disappears due to the ordering on the r i . The main 
contribution to the multiple integral comes from regions where operators are grouped 
pairwise close together. Then the operator product expansion can be used to reduce 
the multi-point function. With a perturbation in the bulk of the system the following 
reduction relations are needed 

a{ ri )e(r 2 ) ~ aa( ri )r^ e^e^) ~ tefoKd"" (9) 

The coefficient b in the second relation vanishes when the system is invariant under 
duality (e — > —e). Then the expansion should be continued one step further. This 
occurs for the Ising model only in the bulk [31] since the surface case does not possess 
the duality symmetry. 

Let us first suppose that r 01 is the smallest distance entering the correlation 
function in (8). Then the first relation in (9) can be used and the integration over r 1 
carried out to give 

rr2 />27r 

aa(0) / dr ir \- u - x ' / /(^)d^ (10) 
Ji Jo 

Since in the marginal case ui = 2 — x E , the first integral contributes a factor lnr 2 . It 
may be checked that when one first contracts other pairs, like s(r 1 )e(r 2 ), the result is 
logarithmically smaller. It follows that, to the leading logarithmic order 



6G%> = g n aS f / \nr 2 ((a(0)a(R)s(r 2 ) ■ ■ ■ e (r n )» [] Z{v l )d\ l (11) 



where Sf is the angular integral in equation (10). The same process can be iterated 
n times leading to 

5G% - 1 [gaS f lni?]" ((a(0)a(R))) (12) 
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Here the n\ is restored through the successive integrations involving increasing powers 
of a logarithm. Finally the perturbation series can be summed giving 

G <J<J {R,g) = G„„{R,Q)B? aS t (13) 

to logarithmic accuracy. The order parameter local scaling dimension can be read out 
of the correlation function. It is indeed non-universal and reads 

x a {g)=x (J -ga\ f{9)&9 + 0(g 2 ) (14) 
Jo 

The same method can be applied step by step to the calculation of the energy density 
correlation function using the second reduction relation in equation (9). The local 
scaling dimension is then 

x £ (g) = x £ -gb f{6)&6 + 0(g 2 ) (15) 
Jo 

Let us now consider the case of an extended perturbation centered on the surface 
of a semi-infinite system [32]. In the direction perpendicular to the free surface 
translation invariance is lost and the operator product expansion has to be modified 
accordingly. When the first point is the origin, located on the surface, the structure 
constants in (9) only acquire an angular dependence 

<r(0)e(r) ~ a{e)a{Q)r~ Xc e(0)e(r) ~ b(6)e(0)r- x * (16) 

where 9 is the polar angle measured from the surface. When the second point goes to 
the surface too, i.e. when 9 ~ r _1 or tt — r~ , the operator product decays with an 
exponent which is the surface scaling dimension x s E = d [33]. Then one expects quite 
generally a{6) ~ b{9) ~ (sin 9f- x - . 

For the Ising model, with x s = 1, one indeed obtains f [32] 

a{9) = - sin 9 b{9) = - sin 6 cos 2 9 (17) 

7T TT 

The bulk values in (9) are a — (27r) _1 and 6 = due to the duality symmetry in the 
2D bulk Ising system [31]. In this latter case the expansion is in powers of g 2 . 

The calculation of the two-point correlation functions then proceeds as above in 
the bulk. The main change is introduced by the angular dependence of the structure 
constants and one gets the following results for the first-order corrections to the local 
exponents 

<(<?) =x' a -g [ V a(9)f(9)d9 + 0{g 2 ) x s e {g) =x s £ -g f b(9)f(9)d9 + 0(g 2 ) (18) 
Jo Jo 

provided the angular integrals are not singular. The same restriction applies to the 
bulk results in equations (14) and (15). 



f The polar angle was measured from the perpendicular direction in [32] 
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4. Conformal aspects 

In this section we assume that conformal techniques can still be used in the presence of 
a general extended marginal perturbation. The conformal mapping of the perturbed 
system, either infinite or semi-infinite, onto a strip will be used in the next section to 
get the perturbed gaps, allowing a comparison with the perturbation results for the 
local exponents. 

Plane- to -cylinder transformation 

We use the mapping w = (L/2ir)lnz of the full plane z — rexp(i#) onto a cylinder 
w = u + iv, — oo < u < +00, < v < L [34]. Under this transformation the marginal 
perturbation A(z) with scaling dimension is changed into [18] 

A(w) = b(z) v *A(z) (19) 

with a dilatation factor given by 



b(z) 



dw 



-1 



dz 

The shape function transforms in the same way since g is invariant. Using 

2nu „ 2ttv . . 

r = cxp — — 6 = — — (21) 
±j Li 

one obtains 

. /27rV*,/27ro 
,v ) = a — f -r- 



A(ti,u) = g j (22) 

where translation invariance along the strip, which is directly linked to the marginal 
behaviour, is essential for the conformal properties. The perturbation is generally 
inhomogeneous in the transverse direction although translation invariance is preserved 
in the case of a radial perturbation [26, 27]. 

4-2. Special conformal transformation 

The infinitesimal special conformal transformation also plays an important role. When 
a system is invariant under this transformation, differential equations for correlation 
functions or profiles can be obtained, leading to the asymptotic behaviour of the first 
and completely determining the second [1, 35]. 
Let the transformation be written as 

z' = z + ez 2 w = z' = r' exp(ifl') (23) 

Up to 0(e), one obtains 

r' = r + er 2 cos 6 6' = 9 + er sin 9 b(r, 9) = 1 - 2ercosfl (24) 

Together with (19), this leads to the following transformed perturbation 

dm.f 



A(r',6') = (1 -2er cos (9)^3— 



1 — er sin 9 ^y^ cot 9 + 



d9 
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The perturbation is shape-invariant when the coefficient of e vanishes, i. e. when 

A(r,g)= . . g „. (26) 
v ' |rsin0|f* v ; 

This corresponds to an extended perturbation decaying either from the surface of 
a semi-infinite system [11] or from a line defect in the bulk [22-25, 21]. Then 
the correlation functions satisfy the same differential equation as in the surface 
case [1]. Their asymptotic behaviour can also be obtained, with defect exponents 
replacing surface ones. Transforming the correlation functions via the plane-to- 
cylinder mapping, the gap-exponent relation and the tower-like structure of the 
spectrum immediately follows, in agreement with exact results for these perturbations 
[18, 21]. 



5. Gap exponent relations 

5.1. Hamiltonian limit and rescaling 

We now consider a temperature-like extended marginal perturbation in the 2D Ising 
model at the bulk critical point. On the cylinder the original perturbation is 
transformed into (22) with <p — e and y — 1. In the case of a semi-infinite original 
system (0 < 9 < ir), the transformed system is a strip with free boundary conditions 
(FBC) at v — and v — L/2 on the cylinder whereas for an infinite original system, 
periodic boundary conditions (PBC) have to be taken at v = L. 

We work on a square lattice in the extreme anisotropic limit, with unperturbed 
two-spin interactions K 1 3> 1 in the time direction along the cylinder axis and K 2 <C 1 
in the transverse one. At the bulk critical point K 2 = K* where K\ is a dual coupling 
satisfying tanh(i^) = cxp(— 2K 1 ). The perturbation is assumed to act only on the 
transverse interactions K 2 . 

In the Hamiltonian limit [36, 7] the system becomes anisotropic with a correlation 
length ratio [37] 

£ 2 _ cosh(2if 2 ) 

Tx ~ cosh^) " 2Kl {27) 

Isotropy can be restored by rescaling the lattice parameter a 1 in the time direction 
to a l — 2K\ measured in units of a 2 in the transverse direction. The row-to-row 
transfer operator may be written as 

T = exp(-a 1 W) (28) 

The Ising Hamiltonian then takes the form 



*~5 



where the cr's are Pauli spin operators. The first part in (29) corresponds to a 
homogeneous system at the bulk critical point, properly normalized to give critical 
excitations with velocity v s — 1 [38] . The perturbation term keeps the same amplitude 
as in the continuum expression (22) since, lengths being measured in units of a 2 , 
exp(— Ti) operates a transfer by one unit length in the time direction and each 
perturbed bond in H is associated with one surface unit on the isotropic rescaled 
system. 
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5.2. Diagonalization of the unperturbed Hamiltonian 

The unperturbed part of the Hamiltonian in (29) is rewritten as a quadratic form in 
fermion operators through a Jordan- Wigner transformation [39] 

N N-l 
n c(P) = - E( C " C « - o) - o E ( C « - C nM+l + + ^(4 - C N ){4 + Cl ) 



2 7 2 

n=l 



(30) 



where the chain length N is assumed to be even in the following. For a semi-infinite 
original system one has to take N = L/2 and P = in (30). If the original system 
covers the whole plane then N = L and P — ±1 is an eigenvalue of the parity operator 
V = cxp(i7r^ n c+c n ) which commutes with H c . 
The Hamiltonian is put under diagonal form 

W c (P)=E £ fc(% + %-l/2) (31) 

k 

through a canonical transformation [40, 41]. The squares of the excitation energies 
for the diagonal fermions in (31), e k = 2|sinfc/2|, can be obtained as the solutions 
of two equivalent eigenvalue problems with normalized eigenvectors, <I> fc and ^ k . The 
wave-vector quantization depends on the boundary conditions, i. e. on P : 

fc(P = +l) = (2p+l)£ fc(P = -l) = 2p£ _^< P <^-1 (PBC) (32a) 



fc(P = 0) = (2p+l)^ TT Q<p<N-l (FBC) (32b) 
With periodic boundary conditions, the normalized eigenvectors are given by 



<& fe (n) = (-iyV-sinfcn * fc (n) = (-1)"^ - cos k ( n + - ) 0<k<n (33a) 



(33b) 



= (-l)"/| cos fcn * fc (n) = (-l) n ^sinfe(n+^ - ^ < fc < 



$ (n) = -* (n) = ^£ = -*_» = J= (P = -1) (33c) 
whereas one obtains 

= (- 1 ) W ^ + I cosfc ( n - **W = ( ~ 1) " +1 v W + T sinfcn (34) 
with free boundary conditions. 
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5.3. Perturbation theory 

Up to first order in the perturbation amplitude, the levels of H c in (29) are shifted by 



(35) 



The relevant levels in the following are the ground state |0) and the states \a) and \e) 
of the unperturbed Hamiltonian which are involved in the calculation of the lowest 
gaps. These states are the lowest ones with nonvanishing matrix elements (0|cr^|cr) 
and (0|cr^<7^ +1 |£) with the ground state. 

With periodic boundary conditions and N even, |0) is the vacuum of H c (+1). 
\e) = 77^77^10), which is even, also belongs to the spectrum of W c (+1) and contains 
the two lowest excitations corresponding to p — and p = —1 in k(P = +1). 
\ a ) = which is odd, belongs to the spectrum of H c (— 1). It contains a 

single excitation with vanishing wave-vector corresponding to p = in k{P = — 1). 
As a consequence this state is degenerate with the vacuum |0)^_ 1 j of H c (— 1). 

In the case of free boundary conditions, the ground-state |0) is the vacuum of 
H c (0), \e) is defined as above with two excitations corresponding to p = and p = 1 
in k(P = 0) while \a) = t]q\0) only contains the lowest one. 

The matrix elements in (35), like the Hamiltonian, are obtained making use of the 
Jordan- Wigner and canonical transformations and read 

D (n) = (0|« +1 |0) = - £ * fc (n)0> fc (n + 1) (36a) 



D» = M« +1 |a) = D (n) + 2* (n)$ (n + 1) (36b) 



D e (n) = (eKo~ x n+1 \e) = D (n) + 2 [* (n)0> (n + 1) + ^(^(n + 1)] (36c) 

In the case of periodic boundary conditions, the final expressions for the boundary 
terms (n — N, n+ 1 = 1) have to be multiplied by — P as in (30) and the eigenvectors 
$ fe and ^ k are those involved in the diagonalization of H C (P) where P is the parity 
of the states in the matrix elements. Then, using (32-34), one obtains 

D (n) = (iVsin^-)" 1 (P = +1) (37a) 
D » = ^(l + cot^) (P = -l) (37b) 



D E (n) -D (n) = ~ sin ^ (P = +1) (37c) 



while for free boundary conditions 



2 . / 11- \ ...._■) 



D a (n) -D (n) = - sin (—)+0(N-J) (38a) 
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D £ {n) D Q (n) = 1 sin (^) cos 2 (H) + 0(Ar" 2 ) (38b) 

According to the gap-exponent relation [34], one expects the scaling dimensions of 
the operator a = a, e to be given by 

x a (g) = lim [25 a ( ff ) - £ (. 9 )] <(. 9 ) = lim ^ [£ Q ( S ) - E (g)] (39) 

iV— »oo Z7T N— >oo 7T 

which, together with (35) and (36), leads to the first-order changes 

*a(g) - *«(«) - ^^T, f (^r) [DJn) D ° in)] (PBC) (40a) 

<(<?) - <(0) = to^-g^f (tJf) [ fl a(«) - A>(»)1 (FBC) (40b) 

n 

In the continuum limit, using (37-38) for large TV, this transforms into 

^(9) = \--tfjMf{0) + O{ 9 2 ) (41a) 

x £ (.g) = l + 0(.g 2 ) (41b) 

for the perturbed scaling dimensions near the source of the inhomogencity in the plane 
and 

<{9) = \ - - ^ Mf{6) sin0 + 0(g 2 ) (42a) 

x'Jg) = 2-^- [ d6 f{6) sin6» cos 2 9 + 0(g 2 ) (42b) 
n Jo 

in the half-plane, in full agreement with the results of section 3. 
6. Elliptic defects 

As an illustration of the the perturbation results one may consider an extended defect 
with elliptic symmetry in the 2D Ising model. The angular dependence 

/(<?) = (sin 2 6» + KCos 2 6>)- 1/2 0<k<1 (43) 

interpolates between the line and radial defects corresponding to n — and n = 1, 
respectively. 

In the infinite Ising system, using (41a) 

xA9) = l--K(VT^) + 0(g 2 ) (44) 
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where K is the complete elliptic integral. The radial defect result [26, 27, 42], 
x r7 (g) = 1/8 — g+ 0(g 2 ) is recovered in the limit k — ► 1. The correction term displays 
a logarithmic divergence in the line defect limit, k — > 0, which is linked to the singular 
behaviour of the perturbation at 9 = and ir. Introducing a cut-off, one gets a jump 
of the magnetic exponent at g — and local order at the bulk critical point for g > 
[21, 31, 43]. 

In a semi-infinite system (43) corresponds to a defect with its main axis along the 
surface. Then, using (42) 



1 Ag arcsin VI — k 



0(g 2 



(45a) 



x»{g) = 2--2- 

7T 



arcsin \J\ — n y/n 

(1 - K) 3 /2 



1 - K 



0(g 2 ) 



(45b) 



In the surface defect limit, k = 0, x"„{g) = l/2-2g+0(g 2 ), x s e (g) = 2-4g+0(g 2 ) while 
for the radial defect, k — ► 1, x%{g) = l/2-4g/ir + 0(g 2 ), x s £ (g) = 2- m g /{3ir) + 0{g 2 ) 
in agreement with know exact results [18, 26]. 

When the main axis of the defect is perpendicular to the surface, i. e. with 



one obtains 



f(6) = (cos 2 6» + k sin 2 6>)- 1/2 0<k<1 



in 



+ 0(g 2 ) 



(46) 
(47a) 



a£(ff) = 2- 



8 5 



(1 - nfl- 



Vl^+1 



0(g 2 ) (47b) 



with the same limits as in (45) for the radial defect (n — 1). When k = 0, i. e. 
with a line defect perpendicular to the surface, x s e (g) = 2 — 8g/n + 0(g 2 ) whereas the 
correction to x s a diverges logarithmically. Then one expects, as for a bulk extended 
line defect, a jump of the surface magnetic exponent at g = and local order for 
enhanced couplings. This could be checked using the techniques of [18, 21]. 



7. Conclusion 

Our main result is the extension of the gap-exponent relations to the case of a general 
marginal extended perturbation in the 2D Ising model. Although this result was 
obtained only up to first order in the defect amplitude, one may conjecture that it 
remains true to all orders, like in the exactly solved limiting cases. 

The shape invariance of the perturbation under the special conformal transforma- 
tion, when it decays like a power of the distance to a line, leads on the cylinder to a 
spectrum containing conformal towers. At least in this case further work should allow 
the determination of the spectrum-generating algebra. 

When the perturbation expansion for the exponents contains divergent terms, one 
expects local order at the bulk critical point and a first-order defect transition when 
the couplings arc enhanced as for an internal line defect. 

Finally let us mention that our results could be extended to other two-dimensional 
conformal systems using conformal perturbation theory [6, 7]. 
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